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ABSTRACT 

We present simultaneous dual-frequency radio observations of Cygnus X-3 during a 
phase of low-level activity. We constrain the minimum variability timescale to be 
20min at 43 GHz and 30min at 15 GHz, implying source sizes of 2-4 AU. We detect 
polarized emission at a level of a few per cent at 43 GHz which varies with the total 
intensity. The delay of ~ lOmin between the peaks of the flares at the two frequencies 
is seen to decrease with time, and we find that synchrotron self-absorption and free- free 
absorption by entrained thermal material play a larger role in determining the opacity 
than absorption in the stellar wind of the companion. A shock-in-jet model gives 
a good fit to the lightcurves at all frequencies, demonstrating that this mechanism, 
which has previously been used to explain the brighter, longer-lived giant outbursts 
in this source, is also applicable to these low-level flaring events. Assembling the data 
from outbursts spanning over two orders of magnitude in flux density shows evidence 
for a strong correlation between the peak brightness of an event, and the timescale 
and frequency at which this is attained. Brighter flares evolve on longer timescales 
and peak at lower frequencies. Analysis of the fitted model parameters suggests that 
brighter outbursts are due to shocks forming further downstream in the jet, with an 
increased electron normalisation and magnetic field strength both playing a role in 
setting the strength of the outburst. 

Key words: X-rays: binaries - shock waves - radio continuum: stars - stars: Cygnus 
X-3 - ISM: jets and outflows 



1 INTRODUCTION 

Cygnus X-3 is one of the few persistently bright radio- 
emitting X-ray binary systems in our Galaxy. It exhibits 
resolved relativistic jets during its occasional giant out- 
bursts, which have been directly i maged on milliarcsecond 
scales dMioduszewski et all l200ll : iMiller- Jones et all 12004 
iTudose et ai] T2007) with very long baseline interferometry 
(VLBI). While the giant outbursts are the most spectacular 
and well-studied type of radio behaviour, the system remains 
active at a lower level for most of the time. IWaltman et al.l 
1 19941 ) identified three different types of variability; normal 
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quiescent emission (60-140 mjy, variable on a timescale of 
months), periods of frequent small flaring (< Uy), and the 
major (> 10 Jy) or intermediate (> Uy) flares, with their 
associated pre-flare quench periods. These classifications 
were based on dual-frequency data from the Green Bank 
Interfero meter, with 10-min ob servations made 3-5 times 
per day. IWaltman et al.l (|l996l ) deduced from the higher- 
than-expected scatter in the quiescent flux densities that 
irregular flux variations continue to be present in this low- 
level state, particularly at the higher frequencies, and data 
at higher time resolution has shown that the so-called "qui- 
escent" state is in fact variable on much shorter timescales 
(Molnar, Reid & Grindlay 1984, 1985, 1988), of order tens 
of minutes. Far from being quiescent, the source is continu- 
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ously active at a low level, aside from the pre-flare periods 
of truly quenched emission. 

Newell, Garrett & Spencer (1998) observed some minor 
flaring episodes (peak fluxes ~ 300 mjy) with the Very Long 
Baseline Array (VLBA) and claimed evidence for superlu- 
minal expansion and contraction of the source. No other 
high-resolution observations have addressed the nature of 
the source during the low-level, active state in which it 
spends the majority of its tim e. The resolved natu re of the 
source in the observations of iNewell et al.l (|l998l ) demon- 
strated that the jets are still active even outside the major 
flaring events. Direct evidence for jet activity outside flaring 
events is also seen in GRS 1915+105, in which an AU-scale 
nuclear jet was resolved with the VLBA (Dhawan, Mirabel & 
Rodriguez, 2000) . At the resolution of the images, the jet was 
observed to be continuous rather than knotty, with a length 
which varied with observing wavelength as ~ 10A cm AU. The 
lightcurves of the observations confirmed that this was an 
optical depth effect. We see emission from the r = 1 sur- 
face, which is further out at lower frequencies. This then 
results in a time delay between flux variations at different 
frequencies, as any injection of relativistic plasma at the jet 
base takes time to propagate out to the point at which the 
jet is optically thi n at the observing f requency (as seen in 
GRS 1915+105 bv lMirabel et al.lll998l ). This also gives rise 
to the smoothing out of the variability at lower frequencies, 
since the observed emission is a convolution over a larger 
region. Measurements of the time delay between different 
observing frequencies can thus help constrain the size scale 
of the jets. 

Being optically thick, such self-absorbed jets tend to 
show very low levels of linear polarization. ICorbel et al.l 
(2000) detected linear polarization at a level of ~ 2 per cent 
in the jets of GX 339-4 while the source was in the low/hard 
X-ray state (believed to be associated with steady, low-level 
jets). The position angle of the electric field vector remained 
constant over more than two years, indicating a favoured 
axis in the system, believed to be aligned with the axis of 
the compact jets. Transient jets, on the other hand, regu- 
larly show significantly higher levels of polarized em i ssion , 
arising from optically thin ejecta (e.g. j Fender et ail Il999l; 
iHannikainen et aj]|2000l : iBrocksopp et all 20071 ). which may 
vary in amplitude and position angle with time, as different 
components dominate the flux density at different times, or 
if the source moves out from behind a screen of ionised elec- 
trons which rotate the position angle via Faraday rotation. 

The relativistically-moving transient ejecta 
observed during X-ray binary outbursts (e.g . 



19781; iMarscher fc Gear! Il985h to X-ray binary systems. 
Fender et al. I l|2004h integrated the internal shock scenario 
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19951 ; Hiellming fc Rupenlll995l ) have previously been mod- 
elled either a s discrete knots of a diabatically-expanding 
plasma (e.g. [van der Laanl 1 19661 ; iHiellming fc Johnston! 
1 19881 ; lAtovan fc Aharonianl 1 19991 ) , or as internal shocks 
within a pre-e xisting steady flow (e . g., Fe nder, Belloni & 
Gallo, 2004). lAtovan fc Aharonianl ifl999l ) found that a 
discrete plasmon model requires continuous replenishment 
of relativistic particles as well as radiative, adiabatic, 
and energy-dependent escape losses. Kaiser, Sunyaev & 
Spruit (2000) argued that such energy-dependent es- 
cape losses were inconsistent with continuous particle 
replenishment via shock acceleration, and adapte d the 
internal shock model originally derived for AGN jets (|Reesl 



into a model explaining the disc-jet coupling over the entire 
duty cycle of an X-ray binary system. In this scenario, as 
a source moves from a hard to a soft X-ray state, the bulk 
Lorentz factor of the pre-existing steady jet increases, giving 
rise to shocks where the highly-relativistic plasma catches 
up with pre-existing slower-moving material downstream in 
the jets. In this work, we focus on the shock-in-jet model, 
since internal shocks form a natural explanatio n for the out- 
burst s in the current disc-jet coupling picture dFender et al.l 
2004), and since shocks provide a mechanism for the 
continuous replenishment of relativ istic particles found 
to be required in a plasmon model IjAtovan fc Aharonianl 
1999). While detailed model fitting for the plasmon scenario 
is deemed to be beyond the scope of this paper, where 
relevant, we will provide brief, generic outlines of how the 
two classes of model differ. 

In Section [2] we present dual-frequency radio 
lightcurves of Cygnus X-3, taken when the source was in 
its normal low-level active state, in order to constrain the 
polarization of the source and investigate the opacity effects 
in the jet (Section[3]). In Section]?] we fit the multifrequency 
lightcurves with a shock-in-jet model, in order to draw com- 
parisons with the properties of the giant outbursts. 



2 OBSERVATIONS 

Cygnus X-3 was observed for 8 hours on 2002 January 25 
(MJD 52299) with the Very Large Array (VLA) in its most 
extended A-configuration, under program code AR458. The 
VLA was split into two subarrays, and the source was ob- 
served simultaneously at 14.940 GHz in one subarray and 
at 43.340 GHz in the other. Occasional snapshots at 1.425, 
4.860, 8.460 and 22.460 GHz were also carried out in order 
to characterise the overall radio spectrum. Observations at 
each frequency were made in the standard VLA continuum 
mode (dual circular polarization in two 50-MHz bands, with 
full polarization products being recorded) . The primary cali- 
brator used was 3C 286 (J 1331+305), and the secondary cal- 
ibrator was J 2007+404. The ab solute flux d e nsity of 3C 286 
was set using the coefficients of iBaars et al.l (1 19771 ). For the 
majority of the observing run, calibrator observations of du- 
ration ~ 1 min were interspersed with 10-20-min scans on 
Cygnus X-3. For 25 min in the middle of the run however, 
fast-switching was used, slewing back and forth between the 
target and the calibrator with a cycle time of 150 s, spend- 
ing 100 s on the target and 40 s on the calibrator. This ap- 
proach was designed to reduce tropospheric phase variation, 
allowing for diffraction-limited imaging, which on the long 
baselines at the high observing frequencies used might not 
otherwise have been possible. The data were reduced using 
standard proc edures within the AIPS data reduction package 
l|Greisenll2003l ). Primary referenced pointing (pointing up on 
a calibrator at a lower frequency to derive the relevant point- 
ing offsets for higher frequencies where the offsets may be a 
significant fraction of the primary beam) was used at both 
43 and 15 GHz, with the times of the pointing scans shown in 
Fig-H] An integration time of 3.3 s was used throughout the 
observing run (except for the pointing scans, where 10-s in- 
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tegrations were mandatory) to allow for high time- resolution 
data editing. 

The weather was good throughout the observations, 
with clear skies and wind speeds of < 3kms -1 , so the point- 
ing accuracy should not have been significantly affected by 
wind loading of the antennas, and the atmospheric opacity 
should have been fairly stable. In order to accurately deter- 
mine the flux density scale, gain curves measured by NRAO 
staff in 200f November were used to correct for gravitational 
deformation of the antennas with changing elevation. Opac- 
ity corrections were made using weather (WX) tables and 
the fitted opacity curves of lButlerl (|2002h . 

Cygnus X-3 was also observed at f 5 GHz by the Ryle 
Telescope, as part of an ongoing monitoring programme 
l|Poolev fc Fenderll 19971 ). The Ryle observations overlapped 
our own VLA observations for 1.6 h at the beginning of the 
VLA run. A systematic discrepancy in the absolute flux 
scales at the two instruments was found, with the VLA mea- 
surements being ~ 25 per cent brighter, exactly as noted by 
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iLindfors et al.l (|2007m . Scaling the Ryle data points by this 
factor brought the two lightcurves into agreement. The Ryle 
Telescope observes a single linear polarization (Stokes I+Q), 
as opposed to the dual circular polarization measured by the 
VLA, although the results in Section \2 . 2 1 show that this dif- 
ference is not significant for Cygnus X-3, and is unlikely to 
be the explanation for the discrepancy. 



2.1 Variability 

The lightcurves of the observations and the spectral index 
a (where S v oc v~ a ) between 43 and 15 GHz are shown in 
Fig.[T] It was found that the flux density of J 2007+404 ap- 
peared to vary slightly during the observations. This can 
almost certainly be attributed to the pointing drifting off, 
since the trend was for the derived flux density to decline be- 
tween pointing observations, especially at the higher of the 
two frequencies. In order to compensate for this effect, the 
J 2007+404 points were all corrected to a mean flux density 
of 1.53 Jy at 43 GHz and 2.83 Jy at 15 GHz, using a multi- 
plicative scaling factor. The Cygnus X-3 flux densities were 
then corrected by linearly interpolating between the gains 
thus derived for the secondary calibrator. 

The source is in the optically thin decay phase follow- 
ing the flare at ~ 15 h, with smaller flares superimposed on 
this steady decrease in flux density. The amplitude of these 
events is greater at 43 GHz than at 15 GHz, and they peak 
10-15 min earlier at the higher frequency, reminiscent of the 
behaviour observed in GRS 1915+105 (e.g. Poolev fc Fender] 
Il997l : iMirabel et aT1ll998l : iFender et al.ll2002D . 



2.2 Polarization 

With the wide parallactic angle coverage of an 8-hour ob- 
servation, it was possible to solve for the "D-terms" (the 
leakage of right circular polarization signal into left circu- 
lar polarization feeds and vice versa) using the calibrator 
J 2007+404. Since 3C 286 was resolved by our combination 
of frequency and array configuration, we used J 2136+006 
as a position angle calibrator (assuming a position angle of 
26.9° at 43 GHz and 70.0° at 15 GHz, derived from mOnitOr- 
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Figure 1. Top panel shows lightcurves of the Ryle observa- 
tions at 15 GHz (dotted circles) and the VLA observations at 
43 GHz (black) and 15 GHz (grey). Also shown are measurements 
at 22 GHz (hollow triangles), 8.4 GHz (hollow circles), 4.8 GHz 
(hollow diamonds) and 1.4 GHz (hollow squares). Oh IAT corre- 
sponds to MJD 52299.0. Vertical arrows on the time axis denote 
times of referenced pointing observations. Bottom panel shows 
the instantaneous spectral index between 43 and 15 GHz, 



ing data taken the following day by NRAO stafl(j). However, 
the polarization position angle for this source can change by 
several tens of degrees in two weeks, so the absolute polar- 
ization position angles derived (the electric vector position 
angles; EVPAs) should be treated with caution. Neverthe- 
less, the magnitude of polarization, P — \/ Q 2 + U 2 (where 
Q and U are the polarized flux densities in the Stokes Q and 
U parameters) is reliable. 

At 15 GHz, no evidence for linear polarization 
at the source position was seen, down to an r.m.s. 
level of 0.22 mjy beam -1 for a single 20-min scan and 
0.06 mJy beam -1 for the full dataset. For a mean 15-GHz 
flux density of 177.0 mjy, this corresponds to a 5<r upper 
limit of 0.17 per cent on the linear polarization. At 43 GHz, 
however, significant linear polarization was detected, at lev- 
els of up to 2.5 mjy. The level of linear polarization initially 
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Figure 2. The polarization parameters of Cygnus X-3 at 43 GHz. 
Top panel shows the full total intensity 43-GHz VLA lightcurve, 
middle panel shows the total intensity lightcurve binned in half- 
hour-long intervals (open circles, dashed line, right hand axis) 
and the linear polarization lightcurve (filled circles, solid line, left 
hand axis) binned to the same resolution. Bottom panel shows the 
electric vector position angle (filled circles, solid line, left hand 
axis) and fractional polarization in per cent (open circles, dashed 
line, right hand axis) as a function of time. In the lower two 
panels, the two data points for the first time bin denote the two 
distinct peaks seen in the polarization image. 



decreases in line with the total intensity, but then rises dur- 
ing and after the flare at ~ 21 h (Fig. [2}. No circular po- 
larization was detected at either frequency to a 5a limit of 
0.6 mjy beam -1 . 

Linear polarization has occasion ally been detected in 
previous giant flare s of Cyg X-3 (e.g ISeaciuist et alj 1 1974 
iLedden et al.l [19761 ) . up to levels of ~ 20 per cent, with 
the fractional polarization increasing with frequency. How- 
ever, the majority of observations have shown no signif- 
icant polarization, with u p per limits of a few per cent 
(e.g. iBraes fc Milevl Il972l : iHiellming fc Balickl 1 19721 ) . so 
there are rela t ively few datasets available for comparison. 
iTudose et al] (|2007l ) used the European VLBI Network 
(EVN) to resolve polarized emission from Cyg X-3, find- 
ing a maximum fractional polarization of 25 per cent at 
the boundary where an ejected knot was interacting with 
its environment. For the derived ro tation measure (RM) of 
— 1233 rad m " 2 (|Ledden et al.lll976l ). we would expect a ro- 
tation of -0.5 rad at 15 GHz, and -0.06 rad at 43 GHz. This 
RM was derived during giant outbursts, when the emitting 
regions had propagated further from the core than those we 
see during these small flaring events. It is likely therefore, 
that there is some extra Faraday rotation in the inner region, 
causing depolarization at the lower frequency either due to 
averaging within the synthesised beam of the observations or 
along the line of sight. However, since the theoretical maxi- 
mum fractional polarization for an optically thick source is 
~ 12 per cent, as compared to ~ 70 per cent for an optically 
thin source, we would in any case expect a lower degree of 
polarization at 15 GHz, where the optical depth is greater. 
The beam and line-of-sight depolarization then reduce the 
observed degree of polarization still further, below our de- 
tection threshold. 



At 43 GHz, the polarization position angle is stable 
to within ~ 20° except during the flare at ~ 21 h, which 
causes a significant, albeit temporary, rotation of the po- 
larization position angle, while the degree of linear polar- 
ization increases slightly. If this flare is due to the forma- 
tion of internal shocks in the flow (see Section 3}, the re- 
sulting compression of the magnetic field at the shock front 
could be responsible for the rotation of the polarization po- 
sition angle and the increase in polarized intensity. Bear- 
ing in mind the uncertainty in absolute position angle cal- 
ibration, we see the mean EVPA evolve from —58° ± 10° 
at the start of the observations, through 76° ± 4° during 
the flare at 21 h, back to —84° ± 15° at the end of the 
observations (as shown in Fig. [2}. For optically-thin emis- 
sion, the measured EVPA should be perpendicular to the 
magneti c field orientation. Previous high-resolution obser- 
vations llMioduszewski et al.|[200ll ; iMiller- Jones et alJliool ; 
ITudose et al.ll2007l ) have shown a north-south jet axis with 
a position angle close to 0° . Thus a field principally aligned 
along the jet axis should have an EVPA of ~ 90°, while 
compression due to a shock viewed sideways should give 
an EVPA of ~ 0°. While the EVPA at the end of the ob- 
servations is consistent with a field aligned along the jet 
axis, there are several possible explanations for the inter- 
mediate position angles seen earlier. In the absence of pre- 
cession, either there is additional Faraday rotation beyond 
the 3° expected from the previously derived RM, the EVPA 
calibration is inaccurate, or we are seeing emission from a 
superposition of emitting regions, some with a longitudi- 
nal field configuration (upstream of the shock front), and 
some from the compressed fields at the shock front. Alter- 
natively, if we are seeing the shock f ront close to face on , 
as has been proposed for Cygnus X-3 (|Lindfors et al.ll2007l ). 
owing t o the small angle between the jet axis and the lin e 
of sight |Mioduszewski et alJliooil ; IMiller- Jones et alJliooj ) , 
we would not expect to see an increase in the degree of or- 
dering of the field. In the absence of data at different fre- 
quencies, or at higher spatial and time resolution, it is not 
possible to differentiate between these competing explana- 
tions. 

The rise in fractional polarization towards the end of 
the observation could be explained in part by the expan- 
sion of the emitting region following the large flare at the 
start of the observations. As the source becomes more opti- 
cally thin, the theoretical maximum fractional polarization 
increases. Furthermore, the density of the stellar wind of 
the companion star decreases as the ejecta move outwards, 
reducing the effect of Faraday depolarization in the wind. 
Alternatively, a change in the relative dominance of a de- 
polarized c ore and highly-polarize d ejecta, as seen in XTE 
J 1748-288 l|Brocksopp et al.ll2007l ). could be behind the in- 
crease in fractional polarization (consistent with the hypoth- 
esis of multiple emitting regions proposed to explain the ob- 
served EVPAs). Multifrequency data or resolved VLBI ob- 
servations would be required to test these hypotheses, while 
higher sensitivity would allow us to make higher-time resolu- 
tion lightcurves in linear polarization, and more accurately 
probe how well the polarized emission followed the total in- 
tensity. Upcoming facilities such as EVLA and eMERLIN 
will provide the required increases in sensitivity and frac- 
tional bandwidth to make such observations. 
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2.3 Imaging 

Imaging at 43 GHz with the VLA in A-configuration gives 
an angular resolution of order 50 mas. This would be suffi- 
cie nt to resolve a rcseco nd-scale extensions such as those seen 
bv iMarti et all feOOll L if present during our observations. 
However, given the time-variable nature of the source, the 
variability had to be removed prior to imaging the source 
in order to search for any potential faint extensions. Imag- 
ing without removing the variability showed no extension 
to an r.m.s. level of 1.8 mjy beam -1 . An automated proce- 
dure was written using ParselTongue, a Python interface 
to AlPS, to remove the time- variable core by making im- 
ages of 1-min time chunks, and subtracting the fitted core 
component in the itw-plane before recombining and imag- 
ing the core-subtracted data from all time intervals. At 
15 GHz, there was no observable extension down to an rms 
of 0.6 mjy beam -1 . The subtraction did not work as well at 
43 GHz, possibly due to atmospheric distortions shifting the 
source position slightly. The central source was not fully re- 
moved, although no evidence for extension was seen down 
to an rms of 0.2 mjy beam -1 . 

The differ ent proper motion measurements quoted in 
the literature dMioduszewski et al.ll200ll ; IMarti et aill200ll ; 
iMiller- Jones et al.ll2004h all predict that a flare would take 
several days to become sufficiently extended to be resolved 
at 15 GHz by the VLA in its A-configuration. The 400-mJy 
flare at the start of our observations could not therefore 
have been resolved. The Ryle Telescope monitoring dat43 
show no flares with 15-GHz flux densities exceeding 250 mjy 
betw een the 2001 September outburst (|Miller- Jones et all 
2004) and the beginning of our observations. The lack of 
extended emission down to our rms limit of 0.2 mjy then 
constrains the e-folding time for the decay of any low-level 
flares to be < 5.2 d, and that for the 2001 September giant 
flare to be < 12.5 d. 

2.4 Spectra 

The overall radio spectrum was sufficiently well-sampled on 
four occasions over the course of the January 2002 observ- 
ing run to generate broadband spectra (Fig. [3]). The spectra 
are consistent with canonical synchrotron spectra affected 
by absorption at the lower frequencies. The overall flux den- 
sity decreases with time and the spectral turnover moves 
to lower frequencies. This is the behaviour expected for the 
decay of the 400-mJy flare at the start of the observations, 
as the ejecta move outwards from the core and expand. The 
lower-level variability seen in Fig. Q] is a secondary effect 
superposed on the general decrease, indicative of multiple 
emitting components in the source. 



3 OPACITY EFFECTS 

As noted in Section f2.1l the lower-frequency lightcurve is a 
smoothed and delayed version of the emission at the higher 
frequency, a clear indication of the existence of opacity ef- 
fects , as commonly ob served in flaring s ources such as AGN 
fe.g. lAller et al.lll985l ). supernovae (e.g. IWeiler et al.lll986h . 

2 http://www.mrao.cam.ac.uk/~guy 
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Figure 3. Broadband radio spectra for the four epochs when the 
lightcurve was well-sampled. The times (IAT hours relative to 
MJD 52299.0) of the observations are noted in the bottom left of 
each plot. Note the turnover moving to lower frequency and flux 
density with time. 

and gamma-ray bursts (e.g. Sod erberg et al.l l2006h . Since 
the optical depth for both free-free absorption and syn- 
chrotron self-absorption decreases with increasing frequency, 
we probe more compact regions at high frequencies, with 
lower-frequency emission being delayed until the source has 
either expanded or moved out from behind the absorbing 
me dium and the optical depth has decreased to order unity. 

iMolnar et al. ( 1984) observed similar low-level activ- 
ity in Cyg X-3 with small flares which peaked later and 
with smaller amplitudes at lower frequencies. They mea- 
sured a delay of 15 ± 5 min between 22 and 15 GHz, in- 
creasing to 239 ± 7 min between 15 and 1.4 GHz. They also 
claimed evidence for a periodicity in the range 4.8-5.1 h, al- 
though such a period icity has not since been detected (e.g. 
Ijohnston et al.l 1986, who note that in f act the amplit ude of 
any periodic phenomenon in the data of lMolnar et al.l is less 
than 30 mjy) . To search for any periodicity in our data, we 
constructed the power spectra of the lightcurves at the two 
frequencies, which are shown in Fig. 2] 

3.1 Power spectra 

Since the data were not evenly sampled in time, in order 
to construct the true power spectra of the lightcurves, we 
implemented a one-dime nsional version of the CLEAN al- 
gorithm, as described bv lRoberts et all (|l987t ) and used in 
standard data reduction procedures for aperture synthesis 
radio telescopes. We ascertained the significance level of the 
peaks in the CLEAN spectra by performing Monte Carlo 
simulations, randomly stripping out half the data points in 
the time series and calculating the CLEAN spectra using 
the remaining data i|Heslop fc Dekkersl 120021 ) . The 95 per 
cent confidence limit ags was set at the point below which 
lay 95 per cent of the sorted, concatenated CLEAN spec- 
tra from all 1000 Monte Carlo iterations. The final CLEAN 
spectra of the 43 and 15 GHz data are shown in Fig.[4] to- 
gether with their respective significance levels 0-95 of 2.1 mjy 
and 0.67 mjy respectively. 

There is no significant power on periods shorter than 
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Figure 4. Power spectra of 15- (dashed) and 43-GHz (solid) data. 
The vertical dotted line shows the 4.8-hour orbital period of the 
system. The labelled horizontal dotted lines show the 95 per cent 
significance levels, 0:95, for the 15- and 43-GHz data. 
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Figure 5. Cross-correlation between the 15-GHz and the 43-GHz 
lightcurves. Negative delays indicate that the 15-GHz lightcurve 
lags the 43-GHz lightcurve. Lag bins are of length 2 min. Inset 
shows the peak of the cross-correlation function for the three dif- 
ferent sections of data, plotted without error bars for clarity. The 
solid line is for the entire data set, the dashed line is for Section 
1, the dotted line is for Section 2 and the dot-dashed line is for 
Section 3. 



about 20 min at 43 GHz, and about 30 min at 15 GHz. These 
minimum variability timescales provide approximate con- 
straints on the size scales for the t — 1 surface of the 
emitting region at the two frequencies, equal to 2.4 and 
3.6 AU respe ctively (assum ing material to be travelling at 
c). At lOkpc (Dickey 1983), this corresponds toasize of 0.2- 
0.3 mas. This may be compared with the predicted scatterin g 
size of 448(^/lGHz)- 2 09 mas (|Mioduszewski et al.l [200ll ). 
which equates to 0.2 mas at 43 GHz and 1.6 mas at 15 GHz. 
At 15 GHz, we are probing regions smaller than can be ac- 
cessed with VLBI imaging observations. 

The power at 15 GHz is uniformly lower than at 43 GHz, 
owing to the smoothing out of the lower-frequency emis- 
sion due to the larger size of the emitting region. The most 
significant peaks in the power spectra correspond to a pe- 
riod of about 2 h at both frequencies. The lowest-frequency 
peaks correspond to the overall decrease in the mean flux 
density with time, and are not sampled sufficiently well to 
correspond to a believable periodicity. There is no signifi- 
cant peak at the orbital perio d (4.8 h) at e ither frequency, 
consistent with the findings of lOglev et al.l (|200ll ). 

3.2 Delay between the frequencies 

To quantify the time delay between the two frequencies, we 
cross-correlated the lightcurves using the 'locall y-normalised 
discre te correlation function' (LNDCF) of iLehar et al.l 
(1992). This accounts for the uneven sampling of the data 
and the non-stationary nature of the time series, preventing 
the correlation function from being dominated by interpo- 
lations (as is the case for the standard correlation function 
method) and avoiding the introdu ction of sampling artifacts 
(|Moles et al.ll 19861 ; |Petersodfl993T ). 

The derived cross-correlation function is shown in 
Fig. [5] It peaks when the 15-GHz emission lags that at 
43 GHz by 9.5 ± 1.0 min, broadly consistent with the 15 ± 
5 min delay between 22 and 15 GHz found bv lMolnar et al.l 
(1984). Splitting the data into three sections (16.7-18.8 h, 



20.2-22.2 h and 22.7-24.2 h UT, as labelled in Fig.[T}, the lag 
was found to decrease with time, from 13.2 ± 1.0 min in the 
first section to 7.0T1.0 min for the final section. The peak in 
the cross-correlation function also becomes narrower, with 
a diminishing amplitude with time, as shown in Fig. [5] This 
is suggestive of decreasing opacity through the flaring se- 
quence, a trend also seen during the multi-flare out burst 
sequence of 1994 February-March ((Fender et al.| [l997). and 
explained as a decreasing amount of entrained thermal ma- 
terial with time over the course of the outburst. 

In an attempt to better understand the relation between 
the two lightcurves, we now attempt to fit the data with a 
shock-in-jet model. 



4 THE SHOCK-IN-JET MODEL 

Shock-in- jet models have been successfully fitted to the 
lightcurves of the 1994 outburst of GRO J 1655-40 
(|Stevens et all 2003) and to a series of small-scale oscilla- 
tion events in GRS 1915+105 (|Tiirler et al.l |2004 ). It was 
suggested that the same model could be used to fit the 
giant outbursts of the latter source, by simply scaling the 
timescales, pe ak flux densiti e s and frequencies of the indi- 
vidual events. ILindfors et al.l (2007) applied the shock-in-jet 
model to Cygnus X-3, finding a good fit to the data from 
both the 2001 giant outburst and a sequence of intermediate- 
level flares from 1994. To ascertain whether such a shock 
model was appropriate to the smaller-scale events we ob- 
served in 2002, and to derive some of the physical parame- 
ters of the outburst for comparison with the stronger flares 
seen in this system, we fitted the radio lightcurves with a 
shock-in-jet model. 

4.1 Outline of the model 

An analytic model for the evolution of a shock wave propa- 
gating down an a diabatically-expan d ing co nical jet was orig- 
inally derived by iMarscher fc Gear! (|1985l ). This model was 
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generalised by iTiirler et al.l (|2000h . and further refined for 
appli cation to microquas ar systems by ITiirler et all (|2004T ) 
and lLindfors et all (j2007) . According to this model, a shock 
wave propagating down a jet evolves through three differ- 
ent phases, in which the dominant energy loss mechanism 
evolves from inverse Compton losses to synchrotron emis- 
sion, and finally adiabatic expansion. In the inverse Comp- 
ton stage, the peak emission frequency decreases while the 
peak flux density rises. Once the photon energy density be- 
comes equal to the magnetic energy density, synchrotron 
losses take over, and the spectral peak moves to lower fre- 
quency at roughly constant flux density. Finally, during the 
adiabatic phase, both the turnover frequency and flux den- 
sity decrease with time. 

To fit the observed lightcurves with this model, each 
peak in the lightcurve is associated with the evolution of 
a single shock front. The fitted start time of an event, To, 
is the time of onset of the shock, which has been found to 
correspond fairly closely to the zero-sepa ration time of indi- 
vidual VLBI compo nents in AGN jets (|Tiirler et all 1 19991 ; 
iLindfors et al.l 2006). All shocks are assumed to be self- 
similar, and the characteristics of an individual event (here- 
after referred to as the specificities) are set only by a shift 
in peak flux density, frequency and duration with respect 
to the average outburst, to minimize the number of free pa- 
rameters of the model. 

A deta il ed d escription of the model is given by 
ITiirler et all (|2000h . although we have incorporated a num- 
ber of re finements in fitting th e light curves presented in 
Section [2] iBiornsson fc Aslaksenl (|2000h modified the shock 
model to make the rise in flux density in the initial Compton 
st age less abrupt than was predicted by the original model 
of lMarscher fe Gearl (|l985T ). With this modification, it was 
found to be diff icult to differentiate the Compton and syn- 
chrotron stages jLindfors et alj|2007f ) , particularly in the ab- 
sence of high-frequency infrared data, so we omitted the syn- 
chrotron stage from our modelling. We assumed a very sim- 
ple adiabatically-expanding, conical and non-accelerating 
(constant Doppler- factor) jet. We also assumed a homoge- 
neous synchrotron source, and that all the emission origi - 
nates from the approaching jet, since ILindfors et ail (2007) 
found that accounting for the receding jet in this source did 
not improve the fit to the major and intermediate outbursts 
of Cygnus X-3, most likely due to the small angle of the flow 
to the line of sight. 

This model dif fers from a standard plasmon model 
|van der Laan| [l966) principally by accounting for the rise 
phase of the flare in which radiative (Compton or syn- 
chrotron) losses dominate over adiabatic losses. In this ini- 
tial radiative phase, the electrons cannot move far from 
the shock front before losing energy. The thickness of 
the emission region behind the shock front then varies 
close to quadratically with the jet radius R, with the ex- 
act scaling depending on t he parameters of the model 
ijBiornsson fc Aslaksenl [2000'). Thus the volume decreases 
more rapidly on tracing the emission back to the apex of 
the jet, and due to the decreased source volume, the flux 
density at optically thin frequencies (in the infrared) does 
not become excessively high at small radii. While we have 
no simultaneous infrared data from 2002 January 25, and 
cannot mandate the use of a shock model as opposed to 
a plasmon model, application of this shock-in-jet model 



to previous outbursts o f Cygnus X-3 |Lindfors et al.l 120071 ) 
and to GRS 1915+105 (|Turler et al.ll2004 ) has shown that 
it can reproduce both the infrared and radio lightcurves. 
The second main difference is that since the shock-in-jet 
model involves an underlying flow expanding in two dimen- 
sions, the electron energy decreases more slowly with time, 
as E oc R~ 2 ^ 3 rather than the R -1 scaling expected for a 
three-dimensional expansion of a discrete plasmon with no 
turbulence. This in turn leads to a slightly less steep decay 
of the synchrotron self-absorption turnover flux density with 
frequency and time. 

4.2 Fitting the lightcurves 

The radio lightcurves were decomposed into a series of 10 
individual events (each representing the development and 
propagation of a shock within the flow), plus an initial de- 
caying outburst representing residual emission from events 
prior to the start of the observing period. All frequencies 
were fit simultaneously to determine both the profile of the 
average outburst and the logarithmic shifts in flux density 
(A log S) , frequency (A log v) and time (A log t ) of the tran- 
sition from the Compton to the adiabatic stage, which deter- 
mine the specificities of the individual events. It was found 
that in the absence of infrared data to constrain the high- 
frequency emission, there was significant degeneracy along 
the adiabatic evolution track. To discriminate between dif- 
ferent weakly-constrained solutions for the specificities of in- 
dividual outbursts, an extra constraint was placed on the fit 
parameters, minimizing the combined scatter in the speci- 
ficities A log S, Alog^ and Alogt. The model used 51 free 
parameters (seven to specify the shape and evolution of 
the average event profile, four to represent residual decay- 
ing emission at the start of the observations, plus four per 
outburst to determine the characteristics of each of the ten 
individual events), fitted to 700 data points. The best fitting 
model, giving a reduced x 2 of 20.5 for the entire dataset, is 
shown for the different observing frequencies in Fig. [6] The 
converged model parameters are given in Table [T] and the 
derived start times and logarithmic shifts of the different 
events are listed in Table [2] 

While inspection of Fig.[S]shows that the model appears 
to reproduce the overall shape of the lightcurves fairly well at 
15 and 43 GHz, the reduced \ 2 f° r the overall fit is relatively 
high. It is possible that the error bars on the data points may 
have been slightly underestimated. The uncertainties on the 
data are mainly statistical, although some extra systematic 
errors of 1 per cent were added to the Ryle Telescope data 
and the lower- frequency VLA data. However, they do not 
account for the 3-5 per cent uncertainty in the absolute flux 
scale at the VLA, since an overall gain error should affect all 
data points at a given frequency and have no effect on the 
relative shape of the lightcurves, whereas the discrepancies 
between the model and the data are in the fine details of 
the lightcurves rather than their overall normalisations. It 
is more likely that at these high frequencies, atmospheric 
variations on a timescale shorter than the switching time 
between calibrator and target caused extra scatter in the 
observed flux density. 

However, the high reduced x 2 value also reflects the 
fact that the model is a fairly simplistic idealisation, which 
does not attempt to reproduce in detail the exact shapes of 
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Figure 6. Best-fitting model lightcurves at each frequency, with the data overplotted as open circles. The 15-GHz Ryle Telescope flux 
densities have been multiplied by a factor 1.25 to bring them into agreement with the 14.9-GHz VLA data. The different outburst events 
are shown by dotted lines. The frequency of each lightcurve is marked in the top right of each panel. The IAT start times of the different 
outbursts are indicated by arrows at the top of the figure. 



individual events, and so is unlikely to be able to account 
for all features of the observed lightcurves. We assumed a 
conical jet, a strong shock, no reverse shock, and did not 
account for counter-jet emission or jet precession. The slope 
of the electron spectrum s was assumed to be constant with 
time, and the ratio, v p /vh, of the peak frequency to the high- 
frequency break of the optically-thin synchrotron spectrum 
was assumed to be the same for all outbursts (albeit evolving 
with time within an individual event). Real magnetohydro- 
dynamical jets are likely to be far more complicated, with 



phenomena such as transverse shocks, stationary knots, and 
Kelvin-Helmholtz instabilities, which we cannot begin to ac- 
count for in our self-similar analytic model. The fit is clearly 
sufficient to validate the zero-order picture of adiabatically- 
expanding synchrotron plasma, although in the absence of 
high-frequency millimetre or infrared data, discriminating 
between the rise phases of a shock-in-jet model and a plas- 
mon model is more difficult, and deriving the details of any 
shock structure is clearly impossible without detailed high- 
resolution VLBI imaging. 
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Parameter Value 

s 2.11 

b 1.19 

t p (h) 0.47 

Up (GHz) 26.9 

S p (Jy) 0.14 

f b (THz) 2.30 

t f (h) 0.11 



Table 1 . The seven model parameters specifying the shape of the 
average outburst, s is the index of the electron energy spectrum, 
b is the scaling of the magnetic field with distance z along the 
jet (B oc z~ b ), S p is the peak flux density at the transition to 
the adiabatic stage, v p is the frequency at which the flux density 
peaks, and t p is the time of the peak relative to the onset of 
the shock at time To. is the frequency of the high-frequency 
break in the spectrum where it steepens from v~ ( s ~ 1 " 2 to v~ s / 2 , 
defined at the time of the peak, and tf is the time at which the 
spectrum begins to flatten to v~ ( s_1 )/ 2 prior to the start of the 
adiabatic stage. 

While the electron energy decays more slowly with ra- 
dius in the shock-in-jet model, the addition of turbulence 
into the plasma in a plasmon model can slow the decay of 
the flux density as the plasmon expands. We can examine 
the expected synchrotron flux density durin g the adiabati c 
decay phase for the two different models (e.g. Longair 1994), 

S v o,R'kB^ x)i \-^ ! \ (1) 

where ft is the normalisation of the electron spectrum, 
N(E) oc kE~ 3 , and B is the magnetic field strength. For 
the shock-in-jet model, B oc R~ b and k oc #- 2 ( s + 2 )/3. Thus, 
for the shock-in-jet model, 

S v o,R [w - 4a - 3b(s+1)]/(i u-^ /2 . (2) 

From the values in Table 1, this scales as S„ oc R~ 159 . For 
the plasmon model, B oc R~ 2+v , where r\ is a factor de- 
scribing the degree of turbulence in the plasma, typically 
of ord er 1, with r\ = corresp onding to no strong turbu- 
lence (|Landau fc Lifshitz|[l963l ). The first adiabatic invari- 
ant, Pj_/B, is conserved during the expansion, where p± is 
the component of the particle momentum, p = E/c, perpen- 
dicular to the magnetic field. If the particle distribution is 
isotropic throughout the expansion, then p\ — 2p 2 /3, such 
that E oc R ( -~ 2+rt ^ 2 . Assuming the number of particles in 
a plasmon, VN(E) dE = VkE~ s dE, is conserved, we then 
find that k oc ir [2+s+,)(1 - s )/2] , such that 

s„ocir (i) - 2 V- (s ~ 1)/2 . (3) 

For the same value of s as fitted for the shock-in-jet model, 
a value rj — 1.24 would then provide the same scaling dur- 
ing the decay phase. Thus we cannot use the decay phases 
of the outbursts to conclusively discriminate between the 
two classes of model. Higher-time resolution polarization 
data would help to discriminate between the two types of 
flow, since the compression and consequent ordering of the 
magnetic field at the shock front at the onset of an event 
would lead to a higher degree of polarized emission and a 
rotation of the mean position angle of the field, from one 
ali gned parallel to the jet to one aligned perpendicular to 
it l|Hughes et al.lfl989al ). While there are hints of a rotation 



Epoch T A log 5 A log;/ A log* 



1994 


2.95 d 


1.53 


-0.59 


1.58 




5.80 d 


0.92 


-0.65 


1.62 




7.25 d 


0.73 


-0.16 


1.55 




9.06 d 


0.53 


-0.29 


2.18 




11.44 d 


0.66 


-0.94 


1.18 








-V.oZ 


1 1 Q 




14.16 d 


0.92 


-1.28 


1.75 




16.19 d 


1.33 


-0.72 


1.06 




18.53 d 


0.92 


-0.94 


1.94 




20.37 d 


1.65 


-0.87 


1.14 




22.89 d 


0.90 


-1.15 


1.63 




23.72 d 


1.20 


-1.28 


1.97 




30.09 d 


0.50 


-0.66 


1.43 


2001 


2.64 d 


1.98 


-0.76 


1.76 




4.51 d 


1.80 


-1.14 


1.94 




4.64 d 


1.70 


-0.73 


1.78 


2002 


12.38 h 


0.43 


-0.42 


0.50 




13.47 h 


0.74 


0.49 


0.08 




1 ft Q*l U 


u.uo 




-u.ou 




16.59 h 


-0.10 


0.57 


-0.50 




16.88 h 


-0.11 


-0.16 


0.11 




18.31 h 


-0.07 


-0.22 


0.26 




20.71 h 


0.04 


0.16 


-0.17 




21.24 h 


-0.39 


0.03 


-0.05 




22.13 h 


-0.28 


-0.12 


0.20 




23.06 h 


-0.31 


-0.28 


0.16 



Table 2. The fitted start times and logarithmic shifts in flux 
density, frequency and time, for the three outbursts of Cyg X- 
3, quoted relative to the average outburst derived for the 2002 
dataset. 



of the EVPA during the flare at 21 h (Fig. [2]), the fractional 
polarization does not increase significantly, and higher time 
resolution (i.e. higher sensitivity) would be required to verify 
this signature of a propagating shock front. 

Nevertheless, owing to the good overall match to the 
shape of the lightcurves, we will assume that the basic zero- 
order properties of this model are correct, and go on to use 
the fitted model parameters to constrain the differences be- 
tween the 2002 events and the previously-studied outburst 
sequences of 1994 and 2001. 

The shape of the fitted average outburst in the three- 
dimensional parameter space defined by peak flux density, 
frequency and time after onset is shown in Fig. [7] Also shown 
are the locations of the peaks of the individual outbursts, 
calculated from the average outburst values and the loga- 
rithmic shifts in time, frequency and flux density of the in- 
dividual events (given in Table [2}. While there is significant 
scatter in flux density, they are well-aligned in the frequency- 
time plane along the locus of the spectral turnover (Fig. [7J), 
suggesting that the turnover always follows a similar track 
in frequency and time. The scatter in flux density suggests 
that another factor is at work in determining the brightness 
of a given outburst, possibly the normalisation of the elec- 
tron spectrum, k. Such an anticorrelation between frequency 
and timescale has been previously noted in 3C 273 (Tiirler, 
Cour voisier & Paltani 199 9) and Cygnus X-3 (Lindfor s et al.l 
l2007l) . but not in 3C279 jLindfors et al.ll2006n . 
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Figure 8. Peak flux density, frequency and time after onset of the 
shock, for each of the individual outbursts of Cygnus X-3 in 1994 
(filled triangles), 2001 (filled circles) and 2002 (filled squares), 
plotted against one another. The open markers denote the posi- 
tion of the average outburst for each dataset. The dotted lines 
show the least-squares fits for the scaling of the different parame- 
ters with one another; Sp <x ip 1 ' 62 , S p oc I'p -0 ' 93 , 



Figure 7. The evolution of the average outburst in frequency- 
time-flux density space, (a) indicates the path followed by the 
self-absorption turnover, and (b), (c) and (d) show the projec- 
tions onto the frequency-time, frequency-flux density and time- 
flux density axes respectively. The contours in (b) are points of 
equal flux density, spaced by 0.3 dex. In (c), the grey lines show 
individual spectra at times spaced by 0.3 dex. In (d), they show 
lightcurves spaced by 0.3 dex in frequency. The thick black lines 
show the path followed by the self-absorption turnover, and the 
grey dots show the time, flux density and frequency of the peak 
emission of the individual outbursts, each of which is represented 
by a self-similar scaling of the average outburst. The vertical black 
lines at the top of (b) and (c) show the observing frequencies. The 
two outlying points in (c) and (d) are the two earliest flares (start- 
ing at 12.38 and 13.47 h respectively), which are constrained only 
by the Ryle Telescope data and the later low-frequency points. 



4.3 Comparison to previous outbursts 

A similar shock model was used by iLindfors et aD (|2007T l 
to fit the lightcurves of the 1994 February-March and 2001 
September outburst sequences of Cygnus X-3. Both flaring 
sequences were significantly brighter than the events ob- 
served in this paper. In order to make valid comparisons 
between the different outbursts, we re-fitted those datasets 
using the average outburst and model parameters found 
from the 2002 January 25 dataset. This provided consis- 
tent specificities for the individual outbursts of the different 
flaring sequences. The new logarithmic shifts, which gave 
reduced-x 2 values for the 1994 and 2001 datasets of 11.8 
and 3.0 respectively (as compared to values o f 10.2 and 1.8 
from the original fits bv lLindfors et all 120071 ). are given in 
Table[2] We note that the sparse infrared data from the 1994 
outburst sequence were not well fitted by these model pa- 
rameters, since there were no infrared data available in 2002 
to constrain the high-frequency spectral evolution. However, 
since our goal was to compare the radio emission from dif- 
ferent outbursts in a consistent fashion, we used the av- 
erage outburst derived from the 2002 data to fit all three 
flaring sequences, owing to the simplicity of the model (an 



adiabatically -expanding conical jet, fo llowing the original as- 
sumptions of iMarscher fc Gearll 19851 ). 

In Fig. [8] we use the specificities of the individual out- 
bursts (Table [2} to plot the times, frequencies and flux den- 
sities of the peak emission against one another for all three 
datasets. The peak frequency is anticorrelated with the out- 
burst timescale, not just for the 2002 data, but also be- 
tween the different outburst sequences. Events that evolve 
faster peak at higher frequencies (the calculated Spearman 
rank correlation coeffi cient is —0.79, chanc e probability < 
2 x 10~ 6 ), as found bv lLindfors et all (120071 ). The timescale 
of the outburst also shows a positive correlation with the 
peak flux density (correlation coefficient 0.69, chance prob- 
ability < 10 -4 ), and there is an anticorrelation between the 
peak frequency and peak flux density (correlation coefficient 
—0.68, chance probability < 2 x 10 -4 ). Thus the shorter, 
higher-frequency-peaking outbursts also tend to be fainter. 
Fitted scalings of the three parameters S p , Vp and t p with 
one another are given in Table [3] 

The frequency-timescale correlation could potentially 
be explained by travel time arguments. If the time of peak 
emission is related by travel time to the distance from the 
jet base, then events occurring closer to the core peak earlier 
and, since they have expanded less, do so at higher frequen- 
cies. Self-similar expansion implies that these events evolve 
fastest. 



4.4 Physical parameters 

Assuming that the time of peak flux density can be approx- 
imated by travel time down the jet, this occurs at a distance 
z p — jf3Sct p downstream, at which point the jet radius is 
= z p tan <j> (assuming a cone opening angle of (f> — 2° in 
all cases). The source angular size is given by f9 src = 2R p /d, 
whe re d is the so urce distance (taken to be lOkpc for Cyg 
X-3; lDickevlll983l ). Given the angular size 6 SIC , turnover fre- 
quency, Up, and flux density, S p , and the electron index s, 
with the assumption of a homogeneous spherical source with 
a power-law electron energy spectrum, then the magnetic 
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field B, and the normalisation of the electron spectrum, k, 
can be calculat ed at the turnover using the equations of 
IMarscherl ([ 19871 ). as previously done for a sam ple of Galac- 
tic and AGN jets bv lTiirler fc Lindforsl (|2007l ). While for a 
shock-in-jet model the source has a cylindrical slab-like ge- 
ometry rather than being spherical, the depth scales with 
R during the self-similar adiabatic phase, so approximating 
the source as a sphere will be in error by the ratio of the 
slab thickness x to the jet radius J?, a factor c; ~ 0.005fc~ 5 
(as derived for the case of3C273; iMarscher fc Gearll 19851 ). 
where fc = ^/(lOOkms" 1 Mpc -1 ) and H is the Hubble 
constant. For the c urrently-acce pted value of fc = 0.732 
jSpergel et all 120071 ), £ ~ 0.024. IMarscher fc Gear! (|l985l ) 
also approximated the density of the shocked region as be- 
ing constant up to a distance x = £R behind the shock 
front, at which point it would abruptly drop to zero. They 
found that this should not introduce significant errors into 
the proportionalities for their equations. The assumption of 
homogeneity of the emitting region was further justified by 
iBiornsson fc Aslaksenl l|2000l ). and should certainly be suf- 
ficient to derive the scalings between physical p arameters 
which we now present. Adapting the equations of IMarscherl 
|l987h to account for the non-spherical geometry, and as- 
suming we see the e mitting region with the shock face-on 
l|Lindfors et al.ll200Tl l. the magnetic field is given by 



B p = 10" s 6(a) 



the electron normalisation is 



5 c — ' 



5 G, 



(4) 



k p = 10 6 n(a)d" 1 C" 1 6» l 



lt-lfl-(4a+7), -(4a+5) Q 2a+3 



x 8 



-2(a + 2) 2a -3 /r-\ 

v 'erg cm , (5) 



and the optical depth to synchrotron self-absorption at the 
turnover frequency is 

r p (a) = c 2 (a)K p B p 3+2a y 2 u-^ +2a y^e slc dS^+ 2 ^^ (6) 

where the source size # S rc is given in mas, the peak emission 
frequency v p in GHz, the peak flux density S p in Jy and the 
source distance d in kpc. a = (s — l)/2 is the spectral index, 
and n(a) and 02(a) are functions depending logarithmically 
on a, with values of 0.27 and 1.2 x 10 17 respectively for 
a = 0.5. b(a) is a slowly-varying function of a with a value 
of 3.2 for a = 0.5. Away from the peak, the optical depth 
scales with frequency as r oc j/ _ ( 2o! + 5 )/ 2 . 

To meaningfully compare the physical parameters for 
different outbursts, they should be scaled to the same posi- 
tion z along the jet (i.e. to a fixed jet radius R). The mag- 
netic field and electron normalisation scale as power laws 
with jet radius, B oc R~ b and k oc R~ k . b is given in Table[TJ 
and fc — 2(s + 2)/3 for an adiabatic jet flow as assumed here. 
From the peak fluxes, frequencies and times of the individual 
outbursts, the magnetic field, electron normalisation, and 
optical depth to synchrotron self-absorption at a frequency 
of 15 GHz were calculated, scaled to the same jet radius of 
1 AU using the derived values of b and fc, and plotted against 
peak flux density, S p , frequency, v p , and tim e, t p in Fig. 



In all cases we assumed a distance of 10 kpc (Dickey 1983), 
a jet speed (3 — 0.63, an inclination angle of 10.5° to the line 
of sig ht and a jet half-opening angle of 2° (|Miller- Jones et all 
120041 ). W hile we note that slightly di fferent parameters were 





s P 


v P 




B(1AU) 


k(1AU) 


-r(lAU) 


S p 




-0.46 


0.91 


[0.38] 


1.58 


2.40 


v P 


-0.93 




-1.36 


[-0.19] 


-2.23 


-2.64 




0.62 


-0.47 




1.61 


[-0.55] 


2.85 



Table 3. Fitted scalings of the peak turnover flux density S p , 
frequency v p and time t p with one another and the resulting 
predictions, via equations [4] [5] and [6] for the scaling of the mag- 
netic field, electron normalisation, and optical depth at 15 GHz to 
synchrotron self-absorption, normalized to 1 AU, with the three 
parameters S p , u p and t p . The numbers are the power-law index 
to which the parameter in the relevant row is raised to give the 
scaling of the parameter in the selected column, such that, for 



instance, 5 D oc u D 



Correlations which a Spearman rank test 



found bv lMioduszewski et al. l|200ll ). we use a single set of 
parameters here for consistency between outbursts. 



show not to be significant at the 95 per cent level are shown in 
brackets. 



There is clearly significant scatter in the individual 
events (due to the strong dependences of B and k on 5 P , v p 
and #src). However, we can use equations[4][5]and[6]together 
with the derived scalings of S p , v p and t p with one another, 
given in Table[3] to show the expected trends in B, k and r. 
These are plotted as dashed lines in Fig. [5] Spearman rank 
correlations show that the only significant trends are for the 
magnetic field to increase with event duration (correlation 
coefficient 0.64, chance probability 4 x 10 -4 ), the electron 
normalisation to increase with event brightness (correlation 
coefficient 0.47, chance probability 0.015) and decrease with 
peak frequency (correlation coefficient 0.47, chance proba- 
bility 0.034), and for the opacity to increase with brightness 
and duration, and to decrease with frequency. The correla- 
tions which are not found to be significant have the shal- 
lowest power law indices (between —1.0 and 1.0), suggesting 
that they are the most easily masked by the scatter in the 
parameters. The correlation of optical depth with event pa- 
rameters is due to the distance along the jet at which the 
peak occurs. The shorter outbursts peak at higher frequen- 
cies and closer to the core, and have already reached their 
peak flux density at 15 GHz by the time the jet reaches a 
radius of 1 AU, so they are optically thin. The large, bright 
outbursts have not yet expanded sufficiently to become opti- 
cally thin at 15 GHz by the time the jet radius is 1 AU. The 
calculated values for the physical parameters of the average 
outbursts of 1994, 2001 and 2002 are listed in Tableg] With 
the derived magnetic field strengths, the synchrotron loss 
timescales at 43 GHz range from 18 d to 70 y. By the time 
the shock has reached the radio photosphere, its size is not 
limited by radiative losses, having instead been set earlier 
when the magnetic field was higher, since when it has been 
determined by adiabatic expansion. Thus we cannot further 
constrain the size of the emitting region using the radiative 
loss timescale. 

To investigate the decreasing opacity over the c ourse 
of the flare sequence suggested bv lFender et~aT1 (|l997T ) and 
seen in our analysis in Section [3] we performed a Spearman 
rank-order correlation test between the start times of the 
individual outbursts of 2002 January 25 and their peak flux 
densities, frequencies and timescales shown in Table [2] The 
only significant correlation was a trend for the peak flux den- 
sity to decrease with start time, with a Spearman rank coef- 
ficient of —0.83 (< 0.3 per cent probability of a chance corre- 
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Figure 9. Magnetic field (upper panels), electron normalisation 
(middle panels), and 15-GHz optical depth to synchrotron self- 
absorption (lower panels) scaled to a jet radius of 1 AU, and plot- 
ted against peak flux density, frequency and time for each of the 
individual outbursts of Cygnus X-3 in 1994 (filled triangles), 2001 
(filled circles) and 2002 (filled squares), together with the average 
outburst properties for each of the three flaring sequences (open 
symbols). The dashed lines show the predictions for the scalings of 
the different parameters with flux, frequency, and time, according 
to the fitted coefficients shown in Table [3] 



Parameter 


1994 


2001 


2002 


< S P > (Jy) 


1.4 


9.3 


0.14 


< u p > (GHz) 


4.5 


3.6 


26.9 


< t p > (h) 


16.8 


31.9 


0.47 


z P (AU) 


201 


381 


5.7 


R P (AU) 


7.0 


13.3 


0.20 


Bp (G) 


0.25 


0.03 


0.13 


B(1AU) (G) 


2.54 


0.56 


0.02 


k p (erg s— 1 cm — 3 ) 


6.5 x 10~ 2 


1.9 x 10° 


2.1 x 10 3 


k(1AU) (erg s_1 cm -3 ) 


1.4 x 10 1 


2.3 x 10 3 


2.5 x 10 1 


r p (15GHz) 


1.6 x IO" 3 


8.2 x 10" 4 


3.8 x 10- 1 


r(lAU)(15GHz) 


5.6 x 10° 


4.3 x 10 1 


4.3 x 10" 4 



Table 4. A comparison of the fitted parameters of the different 
outbursts of Cygnus X-3 (2002 January 25 from t his paper, and 
2001 September and 1994 February-March from iLindfors et al.l 
120071) , and the physical quantities derived from them. For a given 
flaring sequence, the average time, frequency and flux density of 
the peak in the spectrum are < S p >, < u p > and < t p >. This 
occurs at the transition from the Compton to the adiabatic stage, 
a distance z p downstream in the jet where the jet radius is R p . At 
this point, the magnetic field is B P (S P , u p ,t p ), the electron nor- 
malisation is K p (S'p, Up, t p ), and the optical depth to synchrotron 
self-absorption at 15 GHz is T p (Sp, u p , t p )(15 GHz). Extrapolat- 
ing these last three values back to the point where the jet radius 
is 1 AU, for valid comparison between the flaring sequences, gives 
the magnetic field B(1AU), the electron normalisation k(1AU) 
and the optical depth t(1AU)(15 GHz). 



lation). This trend, due to the decay of the flare at the start 
of the observations, is seen in the individual lightcurves of 
Fig. [T] There was no significant correlation of optical depth 
with start time. If the opacity is indeed decreasing over the 
course of a flaring sequence, that opacity is not generated by 
the synchrotron mechanism. It could be free-free absorption 
due to a decreasing pro portion of entrained thermal mate- 
rial (jFender et al.lll997l ). although such a component was 
not included in our modelling, so its effects cannot be quan- 
tified. We would need simultaneous well-sampled lightcurves 
at lower frequencies to investigate its effects. 



4.5 



The differences between large and small 
outbursts 



VLBI observations dMioduszewski et al. 2001; 

iMiller- Jones et al.1 12004 iTudose et al.l l2007t ) have clearly 
demonstrated that the outbursts of Cygnus X-3 vary in jet 
morphology as well as in strength, showing both one-sided 
and two-sided jets, with the dominant component moving 
to the south in some cases and to the north in others. Such 
differences cannot be explained by varying the magnetic 
field strength or electron normalisation at the onset of the 
shocks. Changes in the accretion flow, or in the immediate 
environment of the binary system, possibly due to changes 
in the strength of the stellar wind, are plausible explana- 
tions for such variations in the jet morphology. However, 
bearing such caveats in mind, we can attempt to use our jet 
model fits to gain some insight into the differences between 
large and small outbursts. Varying the jet orientation and 
environmental parameters will add scatter to the observed 
relations, but we can hope to discern general trends. 

We found that the same shock-in-jet model can ex- 
plain the radio lightcurves of Cygnus X-3 during giant out- 
bursts, intermediate flares, and low-level activity, suggest- 
ing that the th r ee di screte variability classes identified by 
IWaltman et al.l (| 19941 ) are similar phenomena, which differ 
only in scale. This indicates that a jet outflow persists into 
the low-level active state in which Cygnus X-3 spends the 
majority of its time, and that the formation of shocks is re- 
sponsible for the observed events in the radio lightcurves. 
Outbursts which peak earlier do so at higher frequencies 
and with lower peak flux densities. The fitted model param- 
eters strongly suggest that the faster, weaker outbursts oc- 
cur closer to the core, whereas the brighter outbursts occur 
further downstream. For a constant jet speed (as assumed 
in our modelling), this implies that shocks forming further 
downstream take longer to expand by a given factor, giving 
rise to a longer timescale for the emission, and a lower peak 
frequency. 

While Fig. [5] has considerable scatter, as expected from 
the above discussion, the observed trends demonstrate that 
the brighter, lower-frequency peaking outbursts have higher 
electron normalisations. The trends with magnetic field are 
not as clear-cut, although the correlation of longer-duration 
outbursts with magnetic field suggests that the brighter, 
lower-frequency peaking outbursts should also have higher 
magnetic field strengths. It is possible that the shallowness 
of these correlations allow them to be masked by the scat- 
ter in the parameters, making them appear less significant. 
If both the magnetic field and electron normalisation are 
higher for stronger outbursts, it suggests that some signifi- 
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cant fraction of the magnetic field could be generated by the 
electrons themselves. However, the dependences of the ob- 
servable parameters S p , v p and t p on the underlying physical 
parameters B and k are clearly complex, and better mod- 
elling is required to disentangle the relative importance of 
the magnetic field and electron normalisation in determining 
the size of the outbursts. 



4.6 Implications for high-resolution imaging 

While the lightcurves can be well fitted by shocks propagat- 
ing down a jet, it is natural to consider whether these events 
in fact correspond to discrete VLBI components. Hughes, 
Aller & Aller (1989b, 1991) found that a model consist- 
ing of shocks propagating down a jet could explain the ob- 
served evolution of the total intensity, polarized intensity, 
and VLBI morphology o f three extragalact i c radio sources. 

Several studies (e.g. iMutel et al.lll990l ; lAbraham et al.l 

1996) have found that peaks in total intensity lightcurves 
correspond to the ejection of new VLBI components, with 
the flux density of the ejected components roughly cor- 
responding to the increase in total flux density. Sub- 
sequently, the decomposition of ligh t curves into d i screte 
flaring events dValtaoia et al.1 1 19991 ; iTiirler et all 1 19991 ; 
iLindfors et al.ll2006T ) has found a good correspondence be- 
tween VLBI component zero-separation dates and the fitted 
start times of those events. However, the component flux 
densities were less well fitted, probably due to an underesti- 
mate by the model of the contribution of the underlying jet. 
The largest sample studied to date, incorporating 27 blazars, 
found that every resolved VLBI component had a corre- 
sponding flare in the total intensity lightcurve, with zero- 
separation dates which corresponded well to the start dates 
of th e flares inferred from the lightcurves l|Savolainen et all 
2002). But again, while the flux densities of the flares cor- 
related with those of the resolved VLBI components, there 
was significant scatter about this relation. 

Flares in the total flux density are also observed which 
do not ap pear to correspond to resolved ejecta on VLBI 
scales fe.g. IPvatunina et al ] |2006l : lLindfors et al.ll2006T ). The 
favoured explanation (reinforced by evidence for unresolved 
co mponents derived from a closure phase analysis performed 
by ISavolainen et all l2002h is that such events do produce 
ejecta, but that these have faded before propagating out 
to the point where they can be resolved. ISavolainen et al] 
(2002) therefore suggested that every large (> 30 per cent 
increase in flux density) flare in the total intensity lightcurve 
corresponded to the formation of a new shock. This raises 
the question of what sets the decay time of the ejecta. If 
the shock takes a long time to become optically thin at the 
observing frequency, the emission will have faded before it 
can propagate out far enough to be resolved. Similarly, if 
there are significant radiative losses (particularly important 
at higher frequencies), the emission will rapidly fade be- 
low the detection limit. In our model, the location of the 
outburst in the three-dimensional frequency-time-flux den- 
sity space will then determine whether or not a given flare 
fades too rapidly to be resolved at the observing frequency. 
We note however, that the model is an oversimplification, 
and that other events such as instabilities in the flow, mag- 
netic field enhancements, or interaction with the surround- 
ings could cause a rebrightening further downstream, allow- 



ing components that would not otherwise be resolved to be 
detected. While we cannot hope to account for such phenom- 
ena without introducing many extra free parameters while 
attempting to model the observed lightcurves, they could 
help to account for the discrepancies between model predic- 
tions and VLBI observations. 

In the case of Cygnus X-3, VLBI jets have been seen 
out to angular distances of 50-100 mas ( Miodusz ewski et al.l 
l200ll ; iMiller- Jones et al.ll2004l ; iTudose et al.ll2007l ). Since the 
peak frequency v p scales with distance (i.e. time) along the 
jet as Zp 0,47 (Table [3}, the low-level flares seen here, with 
peak frequencies ten times lower than those seen in the gi- 
ant flare of 2001, might be expected to be extended at the 
10-20 mas level, assuming that the speeds of propagation 
were similar in both cases. However, the resolved knots seen 
during giant flares were observed well after the peak of the 
outbursts, while the emission was fading. For an outburst 
peaking at flux densities 100 times lower, it would be chal- 
lenging to detect extension at such low levels. Furthermore, 
travel-time arguments suggest that it would take several 
hours for shocks travelling at c to reach angular distances of 
10-20 mas, for a source distance of lOkpc. Indeed, the con- 
straints on the source sizes derived in Section [3] were 0.2- 
0.3 mas, smaller than the scattering disk and at the limit of 
what can be achieved with current VLBI arrays, even at the 
highest frequencies. We therefore find it unlikely that these 
flaring events could be resolved with current instrumenta- 
tion. 



5 ABSORPTION MECHANISMS 

In Section [3] we quantified the opacity effects which led to 
the delay and smoothing of the 15-GHz lightcurve with re- 
spect to that at 43 GHz. However, the nature of such opacity 
was not explicitly determined. Plausib le absorption mecha - 
nisms are synchrotron self-absorption l|Mirabel et al.|[l998l L 
and free- free abs orption from either the wind of the com - 
panion star (e.g. iFender et al.lll995l ; IWaltman et al.lll996l). 
or fro m entrained thermal electrons in the jet ( Fender et al.l 
Il997h . 

We can quantify the effect of the dense wind of the 
Wolf-Rayet companion star, which causes delays as the 
emitting regions move out from behind the radio photo- 
spheres, located closer to the central binar y at higher fre- 
quenc ies. Using the parameters suggested bv lWaltman et aD 
( 1996) for the wind speed, temperature, ionisation param- 
eters and mass-loss rates would predict the 15 and 43 GHz 
photospheres to lie at ~ 4 x 10 13 and ~ 2 x 10 13 cm re- 
spectively, implying a delay of ~ llmin between the two 
frequencies, assuming propagation at the speed of light. 
The speed of the ejecta ha s been suggeste d to lie between 
0.3c JSchalinski et al.lll995h and 0.81c |Mioduszewski et al.l 
l200ll L corresponding to delays of 14-36 min. With only two 
frequencies, we cannot verify whether the relation holds 
across all frequency bands, but it is consistent with the ob- 
served dual-frequency delay only if the speed of the ejecta is 
at the upper end of the allowed range. If this were the only 
source of opacity, the decreasing delay seen over the course 
of the outburst sequence would suggest variations in either 
the speed of the ejecta, or the location of the photospheres 
between outbursts. Since variations in the wind would take 
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> 12 h to propagate out to the photospheres at a wind speed 
of 2000 km s -1 , we consider the latter possibil ity unlikely. 

The frequency-dependent delays found bv lMolnar et al.l 

(| 1984T ) are inconsistent with the predictions of free-free ab- 
sorption in a wind. The ratio of delays between different 
pairs of frequencies can be compared with that expected for 
emitting knots moving at consta nt speed between photo - 
spheres whose radii scale as u~ -°' 67 (Wr ight fc Barlowll 19751) . 
The measured delay between emission at 2 and 20 cm is 
longer than predicted, suggesting the presence of an addi- 
tional source of opacity affecting the 20 cm emission, even in 
small flares of amplitud es comparable to those presented in 
this work. As found by Wa ltman et al. I d!996h . while wind 
opacity may play a role, a variable injection rate of relativis- 
tic electrons into the jet dominates the observed lightcurves. 
If wind opacity is not a significant contributor to the overall 
delay, the photospheric radii must be smaller than calcu- 
lated, implying that the wind mass loss rate is low e r tha n 
the value of 1O _5 M0 y _1 assumed bv lWaltman et al.li| 19961 ). 

The shock-in-jet model of Section [4] is able to fit the 
lightcurves fairly well at both frequencies, and thus is able 
to account for the observed degree of delay and smooth- 
ing. It takes no account of any free-free opacity, so syn- 
chrotron self-absorption alone is certainly capable of gen- 
crating the observed delay. The same model also reproduces 
the giant outbursts, in which the shocks form further down- 
stream in the jet, beyond the radio photospheres formed 
by the wind of the companion star, when synchrotron self- 
absorp tion or free-free abs orption by entrained thermal ma- 
terial (|Fender et al.l 119971 ) is believed to be responsible for 
the frequency-dependent delays. With only two relatively 
high radio frequencies to determine the delay, it is not pos- 
sible to constrain the absorption mechanism further with 
our data. While the decreasing delay with time is suggestive 
of entrained thermal material, t he decrease is not as marked 
as found bv lFender etafl l|l997l) for larger flares. 



6 CONCLUSIONS 

We have observed an episode of low-level activity in the 
X-ray binary system Cygnus X-3, demonstrating that the 
source remains active in its so-called quiescent state. The 15- 
GHz lightcurve lags the 43-GHz lightcurve by 9.5 ± 1.0 min, 
and the variability is smoothed out at the lower frequency. 
There is no significant power on timescales shorter than 
~ 20 min at 43 GHz, suggesting an approximate size scale 
of ~ 2 AU for the emitting region at this frequency. As ex- 
pected, the source was not seen to be extended at either 
frequency. Linear polarization was detected at levels of 1— 
3.5 per cent at 43 GHz, but not at 15 GHz, suggesting that 
either opacity, beam or line-of-sight Faraday depolarization 
is affecting the lower frequency. The polarized flux density 
varies with the total intensity, with an increase in polarized 
flux and a change in the electric vector position angle dur- 
ing a flare in total intensity, consistent with the formation 
of shocks in the flow aligning the magnetic field and giving 
rise to enhanced emission at the shock front. 

The lightcurves could be well fitted at all frequencies 
with the same shock-in-jet model as has been applied to 
larger outbursts, suggesting that the same phenomena are 
at work in causing these low-level flares as are responsible for 



the giant outbursts. Stronger flares are seen to peak at lower 
frequencies and on longer timescales than fainter outbursts, 
consistent with the formation of shocks further downstream 
in the jet. Brighter outbursts which peak at lower frequen- 
cies tend to have higher electron normalisations, and we find 
evidence that the strength of the magnetic field may also be 
higher in such outbursts. However, the scatter in these re- 
lations demonstrates that the true dependencies of the ob- 
servables on the fundamental physical parameters are sig- 
nificantly more complex than a simple power-law relation. 
An in-depth investigation of a plasmon model is beyond the 
scope of this paper, and on the basis of these data alone we 
cannot definitively rule out a plasmon model in favour of 
a shock-in-jet scenario. However, the similarities we found 
to the larger outbursts with better infrared coverage (a fre- 
quency regime where the predictions of the two models can 
in principle be distinguished) provide qualitative support for 
our application of the shock-in-jet model to explain these 
small-scale flaring events. 

In agreement with previous work, we find that free-free 
absorption in the stellar wind of the companion star is un- 
likely to be responsible for the observed opacity effects and 
the time delay between the two frequencies. The decreas- 
ing delay with time is suggestive of a decreasing fraction of 
entrained thermal material, but with only two frequencies, 
we cannot quantify the relative contributions of synchrotron 
self-absorption and free-free absorption from entrained ther- 
mal material. 
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